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Topological phases have been explored in various fields in physics such as spintron- 
ics, photonics, liquid helium, correlated electron system and cold-atomic system. This 
leads to the recent foundation of emerging materials such as topological band insula- 
tors, topological photonic crystals and topological superconductors/superfluid. Here, 
we propose topological magnonic crystal which provides protected chiral edge modes 
for magnetostatic spin waves. Based on a generic linearized Landau-Lifshitz equation, 
we show that a magnonic crystal with the dipolar interaction usually acquires spin 
wave band with non-zero Chern integer. We argue that such magnonic systems are 
accompanied by the same integer numbers of chiral magnonic edge modes, along which 
the spin wave can propagate in a unidirectional manner without being backward scat- 
tered. Being a robust 'one-way spin-wave guide', the proposed chiral magnonic edge 
mode makes it possible the magnonic Fabry-Perot type interferometer, which can be 
utilized as the spin-wave logic gate. 



Topological phases in condensed matters have been at- 
tracting much attention because of their fascinating phys- 
ical properties. Discoveries of topological band insula- 
torgUSElopen U p emerging research paradigm on spin-orbit 
interaction physics. Relativistic spin-orbit interaction in 
the topological band insulator endows its Bloch electron 
bands with non-trivial global structures, which lead to 
novel surface metallic states^ and topological magneto- 
electric effectP Superconductor analogues of topological 
insulators have exotic edge modes^ or bound statesp^ 
which are composed only of 'real- valued' fermion field 
dubbed as Majorana fermion. Some aspe cts o f these 
bound states are experimentally confirmed] 11 * 12 ! foster- 
ing much prospect of the realization of quantum comput- 
ers.^ Photonics analogue of topologi cal ph ase with chiral 
edge modes are proposed theoreticall y 1 14 1 15 1 and are subse- 
quently designed in actual photonics crystals.^ Unidirec- 
tional propagations of electromagnetic wave along these 
chiral edge modes were experimentally observed, which 
provides these metamaterials with unique photonic func- 
tionality. 

Here, we theoretically propose a magnonics analogue 
of topological phases, which has topologically-protected 
chiral edge mode for spin waves. Spin wave is a col- 
lective propagation of precessional motions of magnetic 
moments in magnets. Depending on its wavelength, spin 
waves are classified into two categories. One is exchange 
spin-wave with the shorter wavelength, whose motion is 
driven by the quantum-mechanical exchange interactions 
('exchange-dominated' region). The other is ma gnet o- 
static spin wave with the longer wavelength) 17 ! 18 * whose 
propagation is caused by the long-range dipolar interac- 
tion ('dipolar' region). Magnonics research investigates 
how these spin waves propagate in the sub-micrometer 
length scale and sub-nanosecond time scaleP^HUl Lik e j n 
other solid-state technologies such as photonics, phonon- 
ics and plasmonics, the main application direction is to 
explore ability of the spin wave to carry and process 
information. Especially, the propagation of spin waves 



in periodically mod ulated m agnetic materials dubbed as 
magnonic crystal^lHSEsGIl are of one of its central con- 
cern. Owing to the periodic structuring, the spin wave 
spectrum in magnonic crystal acquires allowed frequency 
bands of spin wave states and forbidden-frequency bands 
(band gap). 

We introduce a topological Chern number in these 
spin-wave bands with a band gap. Based on a simple 
magnonic crystal with wide varieties of material parame- 
ters, we predict that the first Chern number for the low- 
est spin- wave band usually takes non-zero integer values, 
when the wavelength is in the dipolar region. The rele- 
vant length scale turns out to be from sub-/xm to hun- 
dreds of /jm, well within the range of nanosize fabrica- 
tion. The nonzero Chern integer in magnonic crystals 
results in the same integer numbers of topological chiral 
edge modes, along which the spin wave can propagate in a 
unidirectional way. Thanks to the topological protection, 
this propagation is generically free from any static back- 
ward scatterings. The unidirectional motion can be ex- 
perimentally measurable especially in yttrium iron garnet 
(YIG), where the coherence length of magnons is on the 
order of centimeters.^ We argue that these edge modes 
can be easily channelized, twisted, split and manipulated, 
which enables to construct novel magnonic devices such 
as spin-wave logic gate. 

Chern number in boson systems 

To introduce topological Chern number for the magne- 
tostatic spin wave, let us consider a general quadratic 
Bogoliubov-de Gennes (BdG) Hamiltonian for boson 
fields in the two-dimensional momentum space (fc-space); 

with f3^ k = [/3i,fc, • • • , pN,h\- N denotes the number of 
bosonic bands, and the 2N by 2N Hermite matrix ffj. 
is composed of the N by N normal parts (particle-hole 
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FIG. 1: Magnonic crystal with chiral edge modes. Pe- 
riodic array of holes is introduced into YIG, forming a square 
lattice. Iron (Fe) is filled inside every hole. 



channel) and TV by TV anomalous parts (particle-particle 
channel) . 

This bosonic BdG Hamiltonian is diagonalized in terms 
of a para-unitary matrix T k instead of a unitary matrix)221 



T k H k T k 



Eh 



(2) 



with [7L7_ fc ]T. 



[/3t,f3-k], and E k is a diagonal ma- 



trix. From the commutation relation for the 7-field, the 
orthogonality and completeness of the new basis are de- 
rived as 



Tl<r 3 T k = <r 3 , T k <T 3 Tl = a 3 



(3) 



respectively where the diagonal matrix cr 3 takes ±1 in 
the particle/hole space, i.e. [<r 3 ]j m = 5j m Oj with cjj = 
+1 for j = 1, • • • , TV and a - = -1 for j = 7V + 1, • • • , 27V. 
In terms of the para-unitary matrix, the field strength 
(Berry's curvature) Bj and the gauge connection (gauge 
field) (Aj iX , Aj iV ) are introduced in the momentum space 
for j-th magnonic band as; 



A jtV = iTr[T l( T 3 Tl(T 3 (d ku T k )], 



(4) 
(5) 



respectively with j = 1, • • • , 27V. Tj is a 27V by 27V diag- 
onal matrix taking +1 for the j-th. diagonal component 
and zero otherwise. The Chern number (the first Chern 
number) associated with the j-th band is given by the 
momentum integral of the respective Berry's curvature 
over the two-dimensional first Brillouin zone (BZ), 



Cj 



1 

2^ 



d 2 kBi 



(6) 



BZ 



This quantity can be shown to be an integer Cj = 
n (see Supplementary Information for details), which 
corresponds to a number of topological chiral edge 
modes. Note that, due to the commutation relation 
of the boson fields, the Chern integer defined in equa- 
tions ( 4|5|6 1 generally takes a form different from those 



Chern integers previously denned in electronic^ and pho- 



tonic^ 
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systems. In the following, we show that a 



two-dimensional magnonic crystal (MC) with the dipolar 
interaction usually supports spin-wave bands with non- 
zero Chern integers. 

2-d magnonic crystals and dipolar interaction 

The MC considered is a ferromagnetic system with its 
magnetization and exchange interaction modulated pe- 
riodically in the 2-dimensional (x-y) direction. For sim- 
plicity, we assume that the system is translationally sym- 
metric along the z-direction, whereas the subsequent re- 
sults are expected to be similar when the thickness of the 
system in the z-direction becomes finite. The system is 
composed of two kinds of ferromagnets; iron and YIG. 
The unit cell of the MC is an a x x a y rectangle, inside 
which iron is embedded into a circular region, while the 
remaining region is filled with YIG (Fig.[TJ). The uniform 
magnetic field Hq is applied along the z direction, such 
that the static ferromagnetic moment M s in both regions 
is fully polarized in the z direction. Propagation of the 
transverse moments (m x , m„ ) is de scribed by a linearized 
Landau-Lifshitz equatiorP^2Mlll 



1 dm± 
MMo dt 



= ± 2iM s (V • QV)to± =f 2im± (V • QV) M s 



=p iHom± ± ih±M s 



(7) 



with V = (d x , d y ), m± — m x ± im y and h± — h x ± ih y . 
(h x ,h y ) stands for the transverse component of long- 
ranged magnetic dipolar field h, which is related to 
the ferromagnetic moment rn = (m x ,m yi M s ) via the 
Maxwell equation, i.e. Wxh — c~ x dte z and V-(h+m) = 
0. The former two terms in the right hand side of equa- 
tion ([T]) comes from short-ranged exchange interaction, 
where Q denotes the square of the exchange interac- 
tion length. M s and Q take the values of iron inside 
the circular region, while taking the values of YIG oth- 
erwise. The filling fraction of the circular region with 
respect to the total area of the unit cell is parameter- 
ized by /. We restrict ourselves to the wavelength much 
longer than atomic lattice constants of YIG and iron, 
so that the use of this Landau-Lifshitz equation is jus- 
tified. Since the velocities of relevant spin-wave modes 
are much smaller than the speed of light, we can further 
employ the magneto-static approximation, replacing the 
Maxwell equations by V x h — and V • (h + m) = 0; 



hu = -d u #, Af = cL 



OyTTly , 



(8) 

with v = x, y. This in combination with equation ([7]) 
gives a closed equation of motion (EOM) for the trans- 
verse moments. In terms of normalized transverse fields; 



(r) 



m+(r) /3 f (r)- m - {r) 



(9) 



with a proper commutator [/3(r), P'(t')] = S(r — r'), the 
coupled EOM reduces to an equivalent generalized eigen- 
value problem with a Hermitian matrix H k (see Methods 
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for details); 



dt 



(3 k 



(3 k 



(10) 



where (31 denotes a vector composed of the Fourier series 

of /?^(r) with respect to the reciprocal vector G; (3' k = 

[••• ,(3l(G),-- -Jand^r) = E fe)G /3fe(G)e-*( fe+G >. The 
EOM is diagonalized in terms of a para-unitary matrix 
TJ, introduced in equation ^ as 



dt 



7fc 
7-fe 



Eh 



'■>-k 



Ik 
7-fe 



[Ek]j gives the dispersion of the j-th magnonic band 



U = I,-'-)* while [T fc ] aj (a 



, G, • • • ) stands for 



the periodic part of the j-th magnonic Bloch wavefunc- 
tion. Correspondingly, the Chern integer is defined for 
each magnonic band via equations ([6]), Q and 

The topological Chern integer thus calculated always 
reduces to zero, whenever the system considered is ei- 
ther time-reversal symmetric; Hl k = -Hfc, or mirror- 
symmetric with the mirror plane being perpendicular to 
the xy plane, e.g. H (k ^ ky) = H (kji _ ky) . Owing to the 
dipolar interaction, however, the anomalous part of H k 
acquires a complex-valued phase factor which depends 
on the momentum in the chiral way (sec Methods for de- 
tails). As a result, these two symmetries are generically 
absent in the 2-d MCs and the bosonic Chern integer can 
take a non-zero integer value. 

This situation is quite analogous to what the rela- 
ti vistic spin-orbit interaction does in fer romag netic met- 
a ] j33 | 34 | an( ^ topological band insulators! 1 ! 2 ! 4 ! 5 ! Moreover, 
contrary to the relativistic spin-orbit interaction, the 
strength of the dipolar interaction is an experimentally 
tunable parameter in the MCs. 2 Q Namely, the dimen- 
sional analysis on equation Q indicates that, when the 
characteristic length scale of the MC (the unit cell size 
A = ^Ja x a y ) becomes larger than the typical exchange 
length y/Q, the dipolar interaction prevails over the ex- 
change interaction. 




FIG. 2: Chern-integer phase diagram for the lowest 
magnonic band and band dispersions of several lowest 
magnonic bands, a, Chern-integer phase diagram for the 
lowest magnonic band as a function of the linear dimension 
of the unit cell size (A = yja x a y ) and the aspect ratio of the 
unit cell (r = a y /a x ) with / = ty x 10 -2 . The phases are char- 
acterized by the Chern integer of the lowest magnonic band, 
Ci. b, Band dispersions of the lowest three magnonic bands 
with r = 1, A = 0.35/im and / = ty x 10~ 2 . A magnonic band 
gap appears between the first and the second lowest band, 
whose size is on the order of 10 _3 THz. c, Band dispersions 
of the lowest and second lowest magnonic band at the critical 
point which intervenes the phase II (Ci = 2) and the phase 
IV (Ci = 0). Dirac cones with same chiralities are formed at 
the two inequivalent X points. 



Chiral magnonic bands in MC 

We found that the Chern integer of the lowest magnonic 
band, Ci, is always quantized to be 2 for the longer 
A, while the integer reduces to zero for the shorter A 
(Fig. [2]}). The respective quantization is protected by 
a finite direct band gap between the lowest band and 
the second lowest band (Fig. [2Jd) . In the intermediate 
regime of A, these two bands get closer to each other. 
With the four-fold rotational symmetry (r = a y /a x = 1), 
the gap closes at the two AT-points at a critical value 
of A (~ 0.28/zm), where the two bands form gapless 
Dirac spectra (Fig. Without the four- fold symmetry 
(r 7^ 1), the band touching at one of the two X points 
and that of the other occur at different values of A. These 
band-touchings are Pi(7r, 0, A Cj i) and P2(0, 7T, A c> 2) in the 



3-dimensional parameter space subtended by two crystal 
momenta k x and k y and the unit cell size A (Fig. 3). 

As in the fermionic casep^H a band-touching point in 
the 3-c£ parameter space generally plays a role of the dual 
magnetic monopole (charge). The dual magnetic field 
considered is a vector field defined in the parameter space 
(k x , k y ,X), generalized from equation Q as the 'rotation' 
of three-component gauge field Aj = {Aj X ,Aj y ,Aj \); 

Bj^VxAj, (11) 

with V = (df- x ,df. ,d\). Here the third component of 
the gauge field is newly introduced in the same way 
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FIG. 3: Schematic configuration of dual magnetic 
monopoles in the 3-dimensional parameter space sub- 
tended by the crystal momenta (k x ,k y ) and the lin- 
ear dimension of the unit cell size (A = yja x a y ) for 
r = a y /a x > 1. Small green spheres denote the magnetic 
monopoles, which emit the dual magnetic field (blue arrows). 
We only show the first quadrant of the Brillouin zone. Three 
phases found in Fig. 2a with r > 1 are separated by the 
A = A Cj 2 plane and the A = A Cj i plane. A magnetic monopole 
located on each plane means that the 'phase transition' is 
accompanied by the change of the Chern integer. 

as previously done in ([5]); 

A jtX = iTriTjazTitT^dxTk)}. (12) 

j specifies either one of the two magnonic bands which 
form the band-touching. At the band-touching point, the 
dual magnetic field for the respective bands has a dual 
magnetic charge, whose strength is quantized to be 2n 
times integer (see Supplementary Information). 

We found that dual magnetic charges for the lowest 
band at the band touching points at Pi and P2 are both 
+2tt (Fig. 3); V ■ Bt = 2tt<5(A - \ c ,i)8{k x - ir)5{k y ) + 
2ttS(X — \c,2)8{k x )5(k y — n), where A C) i < A Cj 2 for a v > a x 
(r > 1) and A C) i > A c-2 for a y < a x (r < 1). The Gauss 
theorem suggests that, when the unit cell size A is varied 
across one of these band-touching points (A c ,i or A c ,2), 
the Chern integer for the lowest magnonic band changes 
by unity, e.g. 

Ci|a>a c ! - Cik i>a = ^~ [ dn-B l = 1, 
27r J Si 

where S\ is a small sphere enclosing the magnetic charge 
at Pi in the 3-<i parameter space. This leads to the phase 
diagram in Fig. [2^,, which describes the Chern integer of 
the lowest magnonic band as a function of the unit cell 



size A and the aspect ratio r. There are four phases, 
divided by the curves A = A c ,i and A = A c ,2- 

A dual magnetic charge is a quantized object, so that 
it can neither disappear by itself nor change gradually. 
Upon any small change of material parameters, two mag- 
netic charges can only move around in the 3-d parameter 
space. Moreover, in the presence of the spatial inver- 
sion symmetry = H p., their locations are always 
restricted at the X points. As a result, the global struc- 
ture of the phase diagram depicted in Fig. [2^, widely 
holds true for other combinations of material parame- 
ters. In fact, for r = 1 and / = tt x 10~ 2 , we found 
A c = 0.370/xm for iron (circular region) and YIG (host), 
and A c = 0.372/im for cobalt (circular region) and YIG 
(host). When varying the filling fraction for iron (cir- 
cular region) and YIG (host) with r = 1, we found 
A c = 0.27 Vn for / = 4tt x 1(T 2 , and A c = 0.348/zm 
for / = 9ir x 1CP 2 . From these observations, we expect 
that, even if materials in the circular region is replaced 
simply by the vacuum, the MC systems with larger unit 
cell size still belongs to topological non-trivial phases, 
having nonzero Chern integers. 

Chiral magnonic edge mode in MC 

The chiral phases with non-zero Chern integers have chi- 
ral magnonic edge modes, which are localized at the 
boundary with the phase with zero Chern integer (phase 
IV) or the vacuum. These chiral modes carry spin wave 
in a unidirectional manner, whose dispersions go across 
the band gap between the lowest and the second lowest 
band. As an illustrative example, we consider a bound- 
ary (y axis) between the MC in phase III and MC in the 
phase IV, whose Chern integers for the lowest magnonic 
bands differ by unity (Fig. [4^.) . The existence of a chiral 
magnonic edge mode between these two is shown from a 2 
by 2 Dirac Hamiltonian derived near their phase bound- 
ary (A = A c ,i). The Hamiltonian generally takes a fol- 
lowing form (see Supplementary Information), 

H cS = ujqt + k(x)t 3 - iad x r x - ibd y T 2 . (13) 

Tj denotes the Pauli matrices subtended by the two- 
fold degenerate eigenstates with their eigen- frequency uq, 
which are formed by the lowest magnonic band and sec- 
ond lowest one at Pi. a and b are constant material pa- 
rameters which are positive- valued. The difference of the 
Chern integers (Ci) for the two phases is represented as a 
change of the sign of the Dirac mass term /c(x); we hence 
suppose that positive k(x) for x > is for the phase III 
and and negative k(x) for x < is for the phase IV (see 
Fig. 4a); lim x - y ± 00 k(x) = ±Koo- The Hamiltonian has a 
following eigenstatej 37 * 38 * 

ip k (r) oc e lkv e~i K( - X "> dx ' \ . 

which is localized at the boundary [x = 0). In terms 
of the surface crystal momentum k along the y axis, the 
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FIG. 4: Chiral magnonic edge modes, a, geometry of 
the system, be, dispersions for Dirac Hamiltonian with the 
Poschl- Teller potential k(x) — k x tanh(i/d) for Kood = 
0.9a (b) and K, x d — 2.9a (c). The bulk magnonic bands 
('bulk') have a gap, and one chiral edge state runs across the 
gap. In c there are some nonchiral edge states, whereas in b 
there is not. Nonetheless, the number of chiral edge modes is 
one, which is determined solely from the difference between 
the Chern integers for the two phases. 

corresponding eigen- frequency is given by E = uq + bk. 
This connects the lowest magnonic band lying at E < 
^0 — and the second lowest band at E > uq + (see 
Fig. 4b, c). The mode is chiral, since it propagates along 
the boundary in a unidirectional way, = d^E = b. The 
argument so far can be generalized into other situations. 
The generalization suggests that the phase with Ci = 1 
at r > 1 (Phase III) or r < 1 (Phase I) has a chiral edge 
mode at its boundary with vacuum, whose dispersion 
crosses the direct band gap at the (n,0) or (0,7r)-point 
respectively, while the phase with C\ = 2 (Phase II) has 
both at its boundary with vacuum. 

The chiral magnonic edge modes can be easily twisted 
or split by the change of the size (A) and shape (r) of 
the unit cell, which we demonstrate in Fig. 5a,b. In 
Fig. 5a, the MC in the phase II (r = 1) is connected 
with the other MC in the phase III (r > 1), whose Chern 
integer for the lowest band differ by unity. The bound- 
ary between these two MC systems supports the chiral 



magnonic edge mode which runs across the direct band 
gap at (0, it) point. The existence of such an edge mode 
can be shown from the other 2 by 2 Dirac Hamiltonian 
derived near A = A c .2. This means that the two chiral 
magnonic edge modes propagating along the boundary 
of the the MC in the phase II (r = 1) are spatially di- 
vided into two, where one mode goes along the boundary 
of another MC in the phase III (r > 1), while the other 
goes along the boundary between these two MCs. This 
configuration realizes a 'spin- wave current splitter', an 
alternative to those proposed in other geometries.^ 

Discussion 

By calculating a newly-introduced bosonic Chern inte- 
ger for spin wave bands, we argue that two-dimensional 
normally-magnetized magnonic crystal acquires chiral 
edge modes for magnetostatic wave in the dipolar regime. 
Each mode is localized at the boundary of the sys- 
tem, carrying magnetic energies in a unidirectional way. 
Thanks to the topological protection, spin-wave propa- 
gations along these chiral edge modes are generally ro- 
bust against imperfections of the lattice periodicity and 
boundary roughness; they are free from any types of 
elastic backward scatterings with moderate strength.^ 
This robustness makes it possible a magnonic analogue 
of Fabry-Perot interferometer as discussed below. 

The interferometer is made up of a coupled of chiral 
magnonic edge modes encompassing a single topological 
MC (see Fig. 5c). Parts of the MC are spatially con- 
stricted by the hole inside, so as to play a role of the 
'point contact' between these edge modesP^ A unidirec- 
tional spin-wave propagation is induced in a chiral mode 
via an antenna attached to the boundary. The wave is 
divided into two chiral edge modes at a point contact. 
Two chiral propagations merge into a single chiral propa- 
gation at the other point contact. Depending on a phase 
difference between these two, the superposed wave ex- 
hibits either a destructive interference or a constructive 
one, which is detected as an electric signal from the other 
antenna. Note that a local application of magnetic field 
changes the velocity of either one of the two chiral edge 
modes locally as v% = b — > U/. = b + Ab(r). This mod- 
ification results in a phase shift of the wave, which thus 
changes the interference pattern observed in the output 
signal. With th e use of this local magnetic field as an 
external control}- 19 * 42 ! the interferometer can serve as a 
solid-state based magnonic logic gate. 

Though chiral edge modes are robust against static 
perturbations, magnetic energies excited in the edge 
mode decay into either phonon states or other magnon 
states in the bulk via inelastic scatterings. The associated 
decay time or coherence length depends on specific ma- 
terials, and spin wave propagation along the chiral edge 
mode survives only within this coherence length. Due to 
the absence of conduction electrons, however, spin waves 
in magnetic insulators are generally known to have long 
coherence lengths. For example, the coherence length in 
YIG in the magnetostatic regime becomes on the order of 
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FIG. 5: Examples of magnonic circuit made by chiral 
magnonic edge modes, a, Schematic picture of a 'spin- 
current splitter'. Two-fold degenerate chiral magnonic edge 
modes are spatially separated from each other at the interface 
region between two MCs with different shapes of the unit cell; 
one with a square unit cell (r = 1), while the other with a 
rectangular unit cell (r > 1). b, Combination of several MCs 
with different sizes and shapes of the unit cell channelizes 
chiral magnonic edge modes at our disposal. Three MCs are 
used, whose material parameters are taken as those of the 
phases I, II and III. c, Magnonic analogue of the Fabry-Perot 
interferometer. An electric current in an input antenna in- 
duces a chiral magnetostatic wave in the edge mode (yellow 
arrows). The input wave is split into two at the point con- 
tact (PCI). After travelling along different trajectories, these 
two waves make an interference at the other point contact 
(PC2). The output wave is measured as an electric signal via 
the other antenna (output). The phase of the chiral magne- 
tostatic wave in the respective trajectory can be controlled 
by the local magnetic field applied near the edge modes (PS1 
and PS2). 



centimeter ™ Thus, all the characteristic spin- wave prop- 
agations depicted in Fig. 5a,b,c are experimentally real- 
izable especially in sizable MC systems made out of these 
magnetic insulators. Measuring these propagations in the 
space- and time-resolved manner is also experimentally 
possible in terms of Brillouin light scatterings J^l which is 
by itself remarkable and which surely leads to the devel- 
opment of innovative spintronic device in future. 

Methods 



Derivation of equation ( 10 ) and material parameters 

To introduce the Chern integer for the magnonic bands, we 
reduced the Landau-Lifshitz equation Q into an equivalent 
generalized eigenvalue problem with bosonic quadratic Hamil- 
tonian as in equation ( |10| . To do this, notice first that the 
commutator between the transverse moments depends on the 
spatial coordinate, i.e. [rh x (r),rh y {r')] = M 3 (r)S(r — r'). To 
remove this spatial dependence, we normalize the fields as; 

m +( r ) a U \ — m -{ r ) 



P(r) 



p\r) = 



(14) 



=4ia(V • QX7)a/3 - 4i/3(V • QV)a 2 



v/2M s (r) v/2M s (r) 

which leads to a proper commutator of the Holstein-Primakov 
(HP) boson; \J3(r), ft{r')] = S(r - r'). In terms of the HP 
boson fields, the Landau-Lifshitz equation is also properly 
symmetrized as, 

dj3 
lit 

- iHo/3 - iad+ty, 
A* =8+ (a£ f ) + d- (a/?) , d± = d x ± id y 

with a(r) = yfM a (r)/2 and |-y| /j,o being omitted for clar- 
ity. The static magnetization and exchange interaction 
are spatially modulated with a lattice periodicity; a(r) = 

Eg Q(G) e iG '' 



(15) 
(16) 



Eg "(G) 



and Q(r) 



with the recipro- 
cal vectors G. Thus, it follows from the Bloch theorem that 
m(r) and ty(r) take a form; 

/^) = EEM G ) e " l(fe+G) '"> 

fe G 
fe G 

where the fe-summation is taken over the first Brillouin zone. 
In terms of these Fourier modes, an equivalent generalized 
eigenvalue problem with a quadratic Hamiltonian for the HP 
field is derived as 



. d 






" Pk " 


i — 
dt 




= [ 





<T',H k . 



(3k 



(17) 



H 



pi 

(3-k 



Hk 



(3k 

PL 



,Pi(G),---pu-G) r - 

,P- k (-G),---P-k(G) 



(Ts in the right hand side takes ±1 in the particle/hole 
space, which comes from the commutation relation of bosons, 
[Pk,P h ] = 1 a nd [p'l , Pk] = — 1. The comparison be- 



tween equations ( 15|16 1 and equation (171 dictates that Hk 
thus introduced is given by a following Hermitian matrix 



(H. — Hk 



Hk 



cx-a + B k + Hol 



a 

a H 



Ik a 

Bk + Ho 1 



(18) 
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iron 


cobalt 


YIG 


M 3 [|iA/m] 


1.711 


1.401 


0.140 


VQlM 


21.24 


28.90 


183.57 



TABLE I: Saturated magnetization M a and exchange interac- 
tion length y[Q of iron, cobalt and yttrium iron garnet (YIG). 



with 



[<x} G>G , = a(G - G'), [I k ] G)G , = 5 G , G , e - 2 ' e ^ G \ 
[B h ]c.G'=4 Yl a(G - Gi)Q(Gi - Ga)a(Ga - G?) 

x (fc + Gi) ■ (fc + G 2 ) 

-4 ^ Q(Gi)a(G 2 )a(G - G?' - G x - G 2 ) 

Gi,G2 

x (G-G')-(G-G' -Gi), (19) 



and 



<9 fc (Gf) = (fc + G)„ + ^(fc + G) y 
|fe + G| 



(20) 



After taking the summation over Gi and Gi in equation ( 19 1 
one can further decompose [Bk]G,G' m to three parts, 



+ 4i ]T [Qa(9 M a)] G _ G ,(G-G% 
+ 4 ]T [Q(d M a)(d M a)] G _ G/ (21) 

p=x,y 

with 

[Q« 2 ] G = 4 / QMa 2 (r) e - lG -d 2 r, 
[Qa(<9 M a 
[Q(9 M a)(9^a 



lo - 5 

1 =1 
Jg - s 

, 1 



Q(r)a(r)(9 (J a(r)) e tG r d 2 r, 



S 



Q(r)(dMr))(dMr))e- lG r d 2 r, 



where the 2-ci integrals in the right hand side are taken over 
the unit cell and S denotes the area of the cell. 

Notice that, owing to complex-valued phase factors in the 
anomalous part, J*. 7^ Jj*, both the time-reversal symmetry 
and the mirror symmetries within the 2-d plane are absent 
in equation ( |18[ ). Without periodic modulations, a = q1 
and Q = Ql, these phase factors can be erased by a proper 
gauge transformation, /3£ -> /3j and /3_ fc -> (3-kVTk, so 
that both symmetries are recovered. In the presence of the 
periodic modulations, a 7^ al and Q 7^ Ql, however, this is 
not the case and these two symmetries are generally gone. 

The specific MC system considere d in this paper is com- 
posed of two kinds of ferromagnetsj 2 ^- 27 ^ 32 ^ the respective 
(square root of) magnetization and (square of) exchange in- 
teraction length are specified by (ay, Qj) (j — 1,2). The unit 
cell of the MC is square shaped, inside which one of these two 
ferromagnets (an, Qi) are embedded within a circular region, 
while the remaining region is filled with the other (02, Q2)- If 
a(r) has a discontinuity at the boundary between these two 
regions, the last term in equation ( |21[ ), [Q(d l _ l a)(d ll a)], gener- 
ally diverges, since having the second derivative with respect 
to a spatial coordinate. Physically, such infrared divergences 
are removed by a smooth variation of the magnetization at 
the boundary. For simplicity, we interpolate a(r) as a linear 
function of the radial coordinate measured from the center of 
the circular region; 



ai (\r\<Ro) 
a ( r ) = \ai- $=2L(oa - a 2 ) (R < \r\ < R 1 ) . (22) 
a 2 (Ri < \r\) 



A discontinuity in Q(r) is also removed by the same linear 
interpolation. The results in Fig. 2 were obtained from the 
numerics with (r ,n) = (R /X,Ri/\) = (0.10,0.125), while 
(M s ,o,M 3A ) = (1.752,0.194) [fiA/m] and (VQ^,V0~T) = 
(33, 130) [A]! 2 -^ As for the material parameters of iron (Fe), 
cobalt (Co), and YIG, we used the values from the literature 1 ^ 
(Table I). To obtain the para-unitary matrix Tk which diago- 
nalizes Hk, we employed the method based on the Cholesky 
decomposition!^ 
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Supplementary Information for 'Topological chi- 
ral magnonic edge mode in a magnonic crystal' 

Quantization of the Chern number 
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The Chern number associated with the j-th magnonic 
band is constructed from the projection operator^ 

Cj = % ^ J^H^-P^KPj)^^)}, (S23) 

where Pj denotes the projection operator filtering out the 
j-th band at the momentum k. The integral is over the 
first Brillouin zone (BZ) in the two-dimensional k space. 
Equation (3) suggests that the projection operator takes 
a form, 



Pj = T fc I> 3 I> 3 , 



(S24) 



where Tj is a diagonal matrix taking +1 for the j- 
th diagonal component and zero otherwise. With this 
definition, equation (S24| allows that J^j Pj — 1 and 



P P = 

-*- j -*■ m. 



5 P 

J m 3 ■ 



denned in equation (S23) reduces 



to a surface integral of the rotation of the gauge field 
Aj = (Aj tX , Aj tV ) over the BZ, i.e. equations (4-6). 

To see its quantization,!^ notice first that such a sur- 
face integral reduces to zero, provided that the gauge 
field is defined uniquely and smoothly over the whole 
BZ. When any of | [Tk] a ,j I has a zero on the BZ, however, 
the gauge field cannot be determined uniquely over the 
whole BZ. In this case, we need to decompose the BZ into 
two overlapped regions^ (Hi and H2 with Hi U H% = BZ 
and Hi n if 2 = dH\ = —dH^ = T), in such a way that 
one of pTfe] ,j (say [Tfe]i 7) does not have any zero within 
one region (Hi), while another (say [Tfc]2j) has no zero 
inside the other (-£^2)- In the former region, we then take 
the gauge, (a\u 



(l)v 

k,ji 



such that (l\u 



k,ji 



is al- 



( 2*) 

ways real positive, while take another gauge (\u kj )) m 



(2) 

the other (H2), making (2\ul ■) to be always real posi- 



tive. Provided that the j-th magnonic band considered 
is isolated from the others (E k j =/= E^.m^j for any fe), 
these two wavefunctions are always related to each other 
by a certain U(l) transformation, 



l«B> = Ki>e Wfc 



k,jl 



(S25) 



Now that the gauge of \u k a ^) and the gauge of \u k j) 

are uniquely defined in Hi and H 2 respectively, A^l = 

i[ a 3]jj( u< k]\ cr 39k, / \u^) (a = 1,2) are smooth functions 
inside each of these two regions. The Stokes theorem is 
applied separately, so that equation (6) is calculated as, 



3 2 



^ / dk V fc x Af ] + / dk V fc x A 

^ J H X 1 27r JH 2 



(2) 




i Idk-(A^-Af) = li dk -V ,<->,,■ 



(S26) 

with Vfc = (9fc x , dk )■ Two regions share a boundary 
(r), which forms a closed loop. 6k introduced in equa- 
tion (S25) generally have a 2irn phase winding (n = Z), 
when fc winds along this loop. Accordingly, the right 



hand side is an integer, i.e. Cj = Z. In the actual 
numerical calculations of Cj, we employed an algorithm 
based on the 'manifestly gauge-invariant' description of 
the Chern integer W 

Magnetic Monopole and Dirac Hamiltonian 

When two bosonic (magnonic) bands form a band- 
touching point in the 3-dimensional p-parameter space 
with p = (k x , k y ,X), the dual magnetic fields associated 
with these two bands via equations (5,11,12) generally 
have a quantized source of their divergence at the point 
(p = p c ). To see this, notice first that, away from the 
band-touching point, p ^ p c , the gauge invariant dual 
magnetic fields for these two bands are uniquely deter- 
mined, which are therefore divergence-free by their def- 
inition. At p = p c , however, the projection to each of 
these two bands cannot be defined, which endows the 
respective dual magnetic field with some singular struc- 
ture. To study this structure, one can use the degenerate 
perturbation theory in a generalized eigenvalue problem. 
The eigenvalue problem takes a form 



H p Tp — <t 3 T p 



with p = (k x ,k y ,\) and p = (—k x ,—k y ,X). The diago- 
nal matrix er 3 takes +1 for the particle space while takes 
— 1 in the hole space. E p is a diagonal matrix, whose ele- 
ments give bosonic (magnonic) bands and are physically 
required to be all positive definite. H p is a Hermitian 
matrix obtained from a quadratic bosonic Hamiltonian 
considered. We decompose this into the zero-th order 
part and the perturbation part; 



H P = H + (Hp - H ) =H + V p . (S27) 

with Ho = H p=Pc . Suppose that Ho has two- fold degen- 
erate eigenstates \v,j) (j = 1,2) with its eigen-frequency 

M> 0); 

H Q \uj) = er 3 \uj) uj 0) 

where the states are normalized as (u 3 ;|er 3 |u, m ) = 5j m . 
On introducing the perturbation V p , the degeneracy is 
split into two frequency levels. The eigenstate for the 
respective eigen-frequency is determined on the zero-th 
order of p — p c as; 

(S28) 

where T diagonalizes H with Tq<t 3 T = T cr 3 T^ = cr 3 
and a unitary matrix U p diagonalizes a 2 by 2 Hamilto- 
nian V^ff formed by the two-fold degenerate eigenstates; 



T p = ToUp + 0(\p-p c \), 



(Ui|Vp|Mi) (ui\Vp\u 2 ) 
(u 2 \V p \ui) (u2\Vp\u 2 ) 



(S29) 



Accordingly, near p = p c , the dual magnetic field de- 
fined in equations (5,11,12) is given only by this unitary 
matrix; 



V x Aj, Aj = iTi[T UpVUp], 
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with V = (dk,.: dk v , d\), while equation (S29) reduces 
to some 2 by 2 Dirac-type Hamiltonian. From this, we 
can prove the quantization of the dual magnetic charge 
at the band-touching point exactly in the same way as in 
the fermionic systems where the sign and the strength 
of the magnetic charge is determined only by the 2 by 2 
effective Dirac-type Hamiltonian. With a proper gauge 
transformation and scale transformation, the effective 



Hamiltonians at the band-touching points Pj (j = 1,2) 
take a form; 

U cS = ojqTq + (A - A cj -)t 3 + ap x Ti + bp v r 2 , 

with a > 0, b > 0, (p x ,Py) = (k x — ir,k y ) for j = 1 and 
(p x ,Py) = {k x ,k y — 7r) for j — 2. From this, equation 
(13) is derived by the replacement of — > —id^. 
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